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INVESTIGATIONS  OF  SEA  ICE  ANISOTROPY, 
ELECTROMAGNETIC  PROPERTIES,  STRENGTH,  AND 
UNDER-ICE  CURRENT  ORIENTATION 


Austin  Kovacs  and  Rexford  M.  Morey 


INTRODUCTION 

The  in  situ  electromagnetic  (EM)  properties  of  sea 
ice  are  important  in  developing  instrumentation  for  the 
remote  sensing  of  ice  thickness  and,  indirectly,  current 
orientation  at  the  ice/water  interface.  Studies  using  an 
impulse  radar  sounding  system  were  made  in  the  spring 
in  1976  and  1977  (Kovacs  and  Morey  1978)  and  in 
1978  (Kovacs  and  Morey  1979a)  on  the  sea  ice  in  the 
area  of  Prudhoe  Bay,  Alaska.  In  these  studies  it  was 
found  that  when  the  crystal  structure  at  the  bottom 
of  sea  ice  had  a  horizontal  c-axis  with  a  preferred  azi¬ 
muthal  orientation,  this  alignment  was  for  all  intents 
and  purposes  parallel  with  the  short-term  (minutes) 
under-ice  current  measurements  made  at  the  time. 

It  was  observed  that  this  oriented  structure  behaved 
as  an  effective  polarizer  of  transverse  electromagnetic 
waves.  The  resulting  effect  was  shown  to  reduce  or 
eliminate  the  EM  signal  reflection  from  the  ice  bottom 
when  the  antenna  E-field  was  oriented  perpendicular 
to  the  preferred  c-axis  direction  of  the  ice  crystal  plate¬ 
lets.  But  when  the  E-f  ield  was  oriented  parallel  with 
the  preferred  crystal  c-axis  direction  a  strong  signal 
return  was  recorded. 

It  was  also  found,  in  sea  ice  having  a  highly  ordered 
crystal  structure  in  which  the  crystals  have  a  preferred 
horizontal  c-axis  alignment,  that:  a)  the  frequency 
spectrum  of  the  signal  reflected  from  the  ice  bottom 
varied  in  the  horizontal  plane,  b)  the  frequency  shift 
was  found  to  be  related  to  the  ice  brine  volume,  c)  the 
relative  change  in  bulk  dielectric  constant  versus  azi¬ 
muthal  angle  did  not  correlate  with  the  coefficient  of 
anisotropy  (in  other  words  a  travel-time  anisotropy 
was  not  measured,  only  a  reflection  anisotropy),  and 


d)  the  surface  of  the  sea  ice  had  no  definite  anisotropic 
trend  as  determined  from  EM  reflection  measurements. 

This  paper  extends  the  previous  studies  on  the  EM 
properties  of  sea  ice  as  determined  with  an  impulse 
radar  sounding  system  operating  on  the  ice  surface  and 
from  an  airborne  platform,  i.e.  a  helicopter.  The  field 
study  was  made  in  March  1979,  in  the  area  of  Prudhoe 
Bay,  at  the  sites  shown  in  Figure  1 . 


FIELD  PROGRAM 

During  the  April  1978  field  program,  electromagnetic 
sounding  measurements  of  sea  ice  were  made  using  an 
impulse  radar  sounding  system  which  radiated  a  time 
domain  wavelet  of  about  14  ns  duration.  The  center 
frequency  of  the  spectrum  of  the  wavelet  radiated  from 
the  linearly  polarized  broadband  antenna  was  about 
1 25  MHz,  with  the  -3  db  points  of  the  spectrum  at 
about  75  and  1 50  MHz.  This  antenna  was  used  as  both 
transmitter  and  receiver.  In  March  1979,  measurements 
were  made  using  linearly  polarized  antennas  with  a  cen¬ 
ter  frequency  of  the  spectrum  of  the  radiated  wavelet 
of  about  280  MHz,  with  the  -3  db  points  of  the  spec¬ 
trum  at  about  220  and  350  MHz.  One  antenna  was 
used  for  transmission  and  one  as  a  receiver.  The  radi¬ 
ated  time  domain  wavelet  was  about  6  ns  in  duration. 
Both  antennas  were  enclosed  in  one  housing;  the  dis¬ 
tance  between  the  centers  of  the  dipole  antennas  was 
about  30  cm. 

As  in  previous  studies,  radar  measurements  were 
made  with  the  radar  antenna  resting  on  the  ice  surface 
and  also  with  the  antenna  elevated  approximately  1 .7 
m  above  the  ice  surface.  In  addition,  measurements 


Figure  1.  Map  of  1 977-1 979  Impulse  radar  studies.  Sites  1 ,  2  and  Tig  were  1 977  or  1 978  study  locations.  Site  G  was 
a  study  location  in  1977,  1978  and  1979.  Sites  WD-79,  BP,  Tig  79  and  FW  were  also  study  locations  In  1979. 
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were  nude  with  antennas  mounted  on  the  side  of  a 
NOAA  (National  Oceanic  and  Atmospheric  Adminis¬ 
tration)  helicopter.  These  measurements  were  made 
with  the  helicopter  resting  on  the  ice  surface,  during 
which  time  the  antennas  were  supported  approximately 
30  cm  above  the  surface,  and  also  while  the  helicopter 
was  in  flight  at  an  altitude  of  1 0  to  1 5  m. 

Radar  measurements  from  the  sea  ice  surface  and 
from  the  helicopter  were  made  at  a  site  about  1  km 
north  of  the  Prudhoe  Bay  West  Dock  (WD-79).  At 
three  other  sea  ice  sites  sounding  measurements  were 
made  with  the  radar  mounted  on  the  helicopter. 
Measurements  were  also  made  on  the  ice  surface  of  a 
freshwater  lake.  At  the  WD-79  site  the  measurement 
procedure  consisted  of  marking  a  30°  increment  polar 
grid  on  the  ice  surface,  making  the  radar  measurements 
at  each  grid  azimuth,  and  obtaining  an  ice  core  from 
the  center  of  the  grid.  Ice  core  temperature  was  meas¬ 
ured  to  0.1°C  using  a  thermistor  bridge.  Measure¬ 
ments  were  made  at  the  1-cm  and  5-cm  depths  and  at 
each  10-cm  increment  below  5  cm.  With  the  use  of  a 
conductivity/salinity  bridge,  the  salinity  of  the  ice 
was  later  determined  from  the  meltwater  of  the  top 
1 14  cm  of  the  ice  surface  and  each  10-cm  section  of 
the  ice  core.  Visual  observation  of  the  ice  core  at 
about  the  1  -  and  1  Vi-m  depths  and  at  the  ice  bottom 
revealed  ice  crystal  structure  and  c-axis  alignment. 
Preferred  c-axis  azimuth  orientation  was  determined 
with  the  use  of  a  large  compass.  These  determinations 
are  believed  to  be  accurate  to  within  1 10°.  Sea  ice 
brine  volume  was  calculated  from  the  ice  core  temper¬ 
ature  and  salinity  data. 

Radar  measurements  were  obtained  with  the  radar 
antenna  elevated  1 .7  m  above  the  ice  surface,  on  top 
of  a  wooden  structure.  A  4-m-square  metal  screen 
was  set  on  the  ice  or  snow  surface  under  the  antenna. 
The  radiated  wavelet  reflected  from  the  surface  of  the 
screen  was  recorded  on  magnetic  tape  and  later  used 
to  provide  a  known  reference  reflection.  The  screen 
was  then  removed  and  the  radar  reflections  from  the 
snow  surface  and  ice  bottom  were  recorded  at  each 
30°  increment  on  the  polar  grid.  The  snow  was  re¬ 
moved,  and  the  above  sequence  of  measurements  re¬ 
peated.  The  antenna  was  then  placed  on  the  ice  sur¬ 
face  and  the  reflection  from  the  ice  bottom  recorded, 
again  at  the  same  30°  grid  azimuths. 

Measurements  made  with  the  antennas  mounted  on 
the  helicopter  were  made  with  the  helicopter  resting 
on  the  ice.  The  radar  reflections  from  the  ice  surface 
and  bottom  were  recorded.  Additional  readings  were 
made  at  30°  azimuth  increments.  This  was  followed 
by  a  series  of  slow,  1 0-m  altitude  flights  over  the  same 
site,  during  which  the  radar  reflection  of  the  top  and 


bottom  of  the  ice  was  again  recorded  along  the  same 
compass  headings  as  when  the  helicopter  was  sitting  on 
the  ice.  An  ice  core  obtained  at  each  site  was  used  to 
determine,  by  visual  inspection  and  compass  measure¬ 
ment,  the  preferred  c-axis  azimuth  direction  at  the  ice 
bottom. 

From  x-y  plots  (see  Appendix)  of  the  average  of 
ten  scans,  the  two-way  travel  time  of  the  radar  signal 
in  the  ice  and  the  relative  voltage  amplitude  of  the  sig¬ 
nal  reflected  at  the  snow  and  ice  surface  and  the  ice 
“bottom”  were  obtained.  The  effective  bulk  dielectric 
constant  e  of  the  sea  ice  was  calculated  from 


(D 


where  c  =  free  space  electromagnetic  signal  velocity 
D  =  tape-measured  ice  thickness  minus  5  cm 
t  =  two-way  travel  time. 


Five  centimeters  was  subtracted  from  the  measured 
ice  thickness  because  at  the  impulse  radar  frequency 
used  in  this  study  the  electromagnetic  boundary  at  the 
bottom  of  growing  sea  ice  has  been  found  to  be  about 
this  distance  above  the  ice/water  interface  (Campbell 
and  Orange  1974),  i.e.  at  or  just  above  the  open  den¬ 
dritic  ice  platelet  structure  found  on  the  bottom  of 
growing  sea  ice. 

The  effective  wavelet  velocity  (Ve)  in  the  sea  ice 
was  determined  from: 


ir  =  2 0=  c  (2) 

f  yje 

An  analysis  of  the  frequency  spectra  (see  Appendix) 
and  the  relative  voltage  amplitude  of  the  signal  reflected 
from  the  various  interfaces  provided  additional  informa¬ 
tion  about  the  material  being  sounded  and  the  nature  of 
the  interfaces.  The  voltage  reflection  coefficient  of  the 
ice  surface  (ps)  was  determined  by  the  ratio  of  the  rel¬ 
ative  peak-to-peak  voltage  amplitude  of  the  ice  surface 
and  metal  screen  reflections.  The  apparent  dielectric 
constant  (e'r)  of  the  ice  surface  was  then  determined 
from: 


r  '  I  1+PS  I 


The  effective  voltage  reflection  coefficient  of  the 
ice  bottom  (pb)  was  likewise  determined,  then  cor¬ 
rected  for  the  partial  reflection  of  the  energy  at  the 
air/ice  interface,  using: 


Pb  =  PbO-Ps) 


(4) 
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Table  II.  Sea  ice  properly  and  radar  reflection  analysis  determination  (snow-covered  ice). 


Snow/Ice  surface -Antenna  elevated _ Ice  bottom -Antenne  elevated 


tee 

Reflection 

Diet 

Max. 

Rett. 

aniso- 

Reflection 

Diet. 

Signet 

Max. 

Refl. 

aniso- 

Snow 

thick. 

Amp.  Cool. 

const. 

reft 

tropy 

Amp.  Coet. 

const. 

vel. 

reft. 

tropy 

thick. 

tm) 

Oit.  (V)  taj 

w 

dir. 

IK) 

Dir  (V)  tab) 

w 

(mini) 

dir. 

IK) 

dir. 

(cm) 

WD-79  169  0* 

373 

-0  271 

30 

0® 

0  24 

-0016 

46 

0139 

30 

322 

-0.234 

26 

30 

- 

- 

- 

— 

60 

298 

-0.216 

2.4 

- 

- 

- 

- 

90 

3.76 

-0272 

3.0 

90 

026 

-0.017 

46 

0.137 

120 

3.33 

-0  242 

2.7 

120 

045 

-0031 

45 

0141 

150 

312 

-0232 

26 

150 

0.36 

-0026 

4.6 

0140 

Avg. 

-0.245 

2.7 

Avg. 

-0.023 

46 

0.139 

Table  III.  Site  average  bulk  brine  volume  and  minimum  effective  elec¬ 
tromagnetic  propagation  losses. 


Site 

Avg. 

brine 

vol. 

CU 

Min. 

two-way 

loss 

(db) 

ice/water 

interface 

loss 

(db) 

Spreading 

loss 

(db) 

Min. 

two-way 

A  +  S  losses 
(db) 

Min. 

one-way 

A  +  S  losses 
( db/m ) 

WD-79 

47 

-29  0  120° 

-7 

-6 

-16 

-5 

1 

76 

-46  0  30° 

-7 

-6 

-35 

-11 

3 

66 

-43  9  90° 

-7 

-7 

-29 

-8 

fW 

0 

14 

-3.5 

-7 

-3.5 

-0.67 

where  p'b  is  the  corrected  effective  voltage  reflection 
coefficient.  This  coefficient  includes  losses  due  to  geo¬ 
metric  spreading,  attenuation,  and  interface  effects. 
Asa  check  on  the  measured  reflection  coefficient ps 
for  the  lake  ice,  the  theoretical  voltage  reflection  co¬ 
efficient  at  the  air/ice  interface  (R^)  was  calculated 
from: 


R 


»li  ~ 


VijTe^T 

Vcj/e,  + 1 


(5) 


where  e,  and  e,  are  the  dielectric  constants  for  the  ice 
surface  and  air,  respectively. 

Similarly,  the  theoretical  voltage  reflection  coeffi¬ 
cient  at  the  ice/water  interface  was  determined  from: 


*i/w 


Vew/el-1 

V^/Cj  +  1 


(6) 


where  ew  is  the  dielectric  constant  of  fresh  water. 

Reflection  coefficients  p$  and  pb  were  used  to  cal¬ 
culate  signal  losses  at  the  ice  surface  and  ice  bottom 
in  db  using: 

|p|  (db)  *  20  logio|p|.  (7) 


Finally,  the  center  frequency  of  the  spectrum  of 
the  time-amplitude  wavelet  reflected  from  the  metal 
screen,  the  snow  and  ice  surfaces  and  the  ice  bottom 
was  determined  using  a  digital  signal  analyzer  (see 
Appendix). 


RESULTS  AND  DISCUSSION 

Radar  reflection  results  for  the  Prudhoe  Bay  West 
Dock  sea  ice  site  in  1979  (WD-79)  and  for  the  fresh¬ 
water  lake  ice  site  (FW)  are  listed  in  Tables  I— III,  along 
with  simitar  sea  ice  results  for  sites  1  and  3  studied  in 
1978  (Kovacs  and  Morey  1979a).  In  addition,  the 
WD-79  data  are  presented  on  polar  coordinate  plots 
from  which  a  coefficient  of  anisotropy  ( K )  is  obtained, 
as  determined  by  the  ratio  of  the  major-to-minor  axis 
of  the  polar  plot  of  the  wavelet  reflection  amplitudes 
vs  E-field  azimuth  (see,  for  example,  Figure  5). 

When  the  antenna  was  elevated,  the  reflection  am¬ 
plitude  of  the  signal  from  the  snow  and  ice  surface 
changed  slightly  with  antenna  E-field  azimuth  orienta¬ 
tion.  These  small  variations  indicated  no  near -surface 
preferred  crystal  orientation,  and  therefore  no  aniso¬ 
tropy  (Fig.  2a  and  2b).  This  was  further  verified  by  our 
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Figure  Z  Relative  voltage  amplitude  of  reflected  wavelet  from  the  snow  (a)  and  Ice  (b)  surfaces,  and  the 
center  frequency  of  the  spectrum  of  the  reflected  signal  from  the  snow  (c)  and  ice  (d)  surfaces  when  the 
antenna  was  elevated  at  site  WO-79. 


visual  inspection  of  the  ice  surface  crystal  structure 
and  by  more  detailed  petrographic  studies  of  sea  ice 
by  many  other  investigators,  e.g.  Langhorne  (in  press) 
and  Weeks  and  Gow  (1979).  Similarly,  polar  plots  of 
center  frequency  vs  E-field  azimuth  orientation  do  not 
exhibit  any  frequency  shift  (Fig.  2c  and  2d).  These 
results  agree  with  those  of  Kovacs  and  Morey  (1 978), 
who  found  that  “all  data  gathered  to  date  using  im¬ 
pulse  radar  indicate  that  the  surface  of  sea  ice  is  either 
not  anisotropic  in  the  horizontal  plane  or  only  weakly 
anisotropic." 

The  center  frequency  of  the  spectrum  of  the  re¬ 
flected  signal  from  the  metal  sheet  was  found  to  be 
about  280  MHz.  However,  the  center  frequency  of 
the  reflected  signal  from  die  snow  and  ice  surfaces 
shifted  down  to  about  200  MHz,  as  shown  in  Figures 


2c  and  2d.  Evidently  frequency-dependent  absorption 
effects  reduced  the  amplitude  of  the  higher  frequencies 
in  the  reflected  wavelet.  Here,  too,  the  polar  plots 
show  that  while  the  center  frequency  varied  slightly 
vs  E-field  orientation,  there  was  no  preferred  azimuthal 
trend. 

The  same  is  true  for  the  center  frequency  of  the 
wavelet  reflected  from  the  ice  bottom  when  the  anten¬ 
na  was  elevated  and  when  it  was  resting  on  the  ice 
surface  (Fig.  3a  and  b).  These  results  indicate  that  the 
frequency-dependent  properties  of  sea  ice  are  not  de¬ 
pendent  upon  azimuth  direction.  However,  the  results 
of  Kovacs  and  Morey  (1 979a)  showed  a  frequency  de¬ 
pendence  vs  E-field  azimuth  at  the  lower  frequency 
wavelet  spectrum  used  in  their  study.  We  do  not  folly 
understand  the  differences  between  these  two  results. 
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Figure  3.  C-axis  direction  and  center  frequency  of  the  spectrum  of  the  reflected  signal  from  the  Ice  bottom 
vs  antenna  E-fleld  azimuth  at  site  WD-  79,  when  the  antenna  was  elevated  (a)  and  when  It  was  resting  on  the 
ice  surface  (b). 


Figure  4.  C-axis  direction  and  relative  voltage  amplitude  of  the  wavelet  reflected  from  the  ice  bottom  w 
antenna  E-fleld  azimuth  at  site  WD-  79,  when  the  antenna  was  elevated  and  the  Ice  surface  was  snow- 
covered  fa)  and  snow-free  (b). 


The  significant  difference  between  the  two  plots 
shown  in  Figure  3  is  the  average  of  the  center  frequen¬ 
cies  of  the  signal  reflected  from  the  ice  bottom:  1 74 
MHz  when  the  antenna  was  elevated  and  1 31  MHz 
when  the  antenna  was  resting  on  the  ice  surface.  The 
latter  is  the  result  of  antenna  loading  effects  which 
occur  when  an  antenna  is  brought  in  contact  with 
another  material.  When  this  occurs,  the  beam  pattern 
and  frequency  spectrum  of  the  radiated  signal  are 
modified. 

The  relative  voltage  amplitude  of  the  reflected  sig¬ 
nal  from  the  ice  bottom  vs  E-field  direction  when  the 
antenna  is  elevated  is  shown  in  Figures  4a  and  b.  Also 
rfwwn  is  the  preferred  ice  bottom  crystal  c-axis  azi¬ 
muth  orientation.  The  variation  in  the  relative  voltage 


amplitude  vs  antenna  E-field  orientation  is  striking;  it 
is  clear  that  maximum  reflection  occurs  when  the  E- 
field  is  aligned  parallel  with  the  predominant  c-axis 
direction.  This  is  in  keeping  with  earlier  findings 
(Kovacs  and  Morey  1978, 1979a).  In  light  of  such 
data,  it  was  deduced  that  the  sea  ice  behaved  as  an  ef¬ 
fective  polarizer  of  the  radiated  electromagnetic  energy, 
and  is  thus  anisotropic.  The  lack  of  measurable  signal 
level  at  the  30-210°  and  60-240°  azimuth  angle  pre¬ 
cludes  construction  of  an  ellipse  which  may  have  taken 
the  form  of  that  in  Figure  5.  However,  an  apparent  co¬ 
efficient  of  anisotropy  was  determined  from  the  ellipse 
drawn  through  the  available  data,  as  shown  in  Figures 
4a  and  b.  The  resulting  apparent  coefficient  of  aniso¬ 
tropy  is  nearly  the  same  (Tables  I  and  II,  WD-79, 
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Figure  5.  C-axis  direction  and  relative  voltage  amplitude 
of  the  wavelet  reflected  from  the  Ice  bottom  vs  antenna 
E-fleld azimuth  at  site  WD-79,  when  the  antenna  was 
resting  on  the  ice  surface. 


antenna  elevated),  with  or  without  a  snow  cover  on 
the  ice  surface:  2.0  and  2.4,  respectively.  The  rela¬ 
tive  voltage  amplitude  of  the  wavelet  reflected  from 
the  ice  bottom  was  slightly  lower  when  the  ice  sur¬ 
face  was  covered  with  3’A  cm  of  snow. 

The  relative  voltage  amplitude  of  the  reflected  sig¬ 
nal  vs  antenna  E-field  azimuth  orientation  when  the 
antenna  was  resting  on  the  ice  surface  is  shown  in 
Figure  5.  Again,  it  is  clear  that  maximum  signal  re¬ 
flection  amplitude  occurs  when  the  antenna  E-field 
is  aligned  parallel  with  the  preferred  c-axis  direction. 

The  last  set  of  radar  measurements  at  the  West 
Dock  site  was  made  with  two  antennas  fixed  to  the 
NOAA  helicopter,  as  shown  in  Figure  6a.  At  this  site, 
only  information  from  one  of  the  two  antennas 
shown  mounted  to  the  helicopter  was  recorded  for 
later  analysis.  Measurements  were  made  while  the 
helicopter  was  on  the  ice  surface  (Fig.  6a)  and  in  flight 
(Fig.  6b).  The  relative  voltage  amplitude  results  for 
the  airborne  antenna  vs  those  for  when  the  antenna 
was  resting  on  the  ice  surface  are  shown  in  Figure  7. 
The  agreement  between  the  airborne  measurements 
and  the  static  on-ice  measurements  is  apparent,  demon¬ 
strating  for  the  first  time  that  it  is  possible  to  detect 
the  existence  of  a  preferred  c-axis  orientation  at  the 
“bottom”  of  sea  ice  from  an  airborne  platform.  Be¬ 
cause  the  c-axes  are  believed  to  be  aligned  with  the 
long-term  current  direction  at  the  ice/water  interface 
(Kovacsand  Morey  1979a,  Weeks  and  Gow  1979, 
Langhorne,  in  press),  it  should  be  possible  to  deter¬ 
mine,  indirectly,  this  azimuthal  alignment  from  an 
airborne  platform  in  remote  areas. 

The  difference  in  relative  voltage  amplitude  for 
the  three  polar  plots  shown  in  Figure  7  is  the  result 
of  various  system  gain  settings  used  for  each  of  the 
measurements.  The  coefficients  of  anisotropy  of  the 
reflected  signal  for  die  antenna  resting  on  the  ice  sur¬ 


face,  mounted  on  the  helicopter  with  the  helicopter  on 
the  ice,  and  airborne  are  3.4, 5.2  and  2.5,  respectively. 
These  results  indicate  that  the  sea  ice  anisotropy,  and 
therefore  the  prevailing  current  alignment  at  the  ice/ 
water  interface,  may  be  more  difficult  to  detect  from 
an  airborne  platform  than  from  the  surface. 

Examples  of  the  impulse  radar  signal  data  collected 
during  a  helicopter  flight  and  displayed  on  a  graphic 
recorder  are  shown  in  Figures  8a  and  b.  The  long 
period  variation  in  the  record  is  due  to  gradual  changes 
in  helicopter  altitude.  The  radar  signatures  from  the 
ice  surface  and  bottom  when  the  antenna  E-field  is 
flown  “parallel”  to  the  preferred  c-axis  direction  are 
apparent  (Fig.  8a).  The  reflected  signal  from  the  ice 
bottom  when  the  antenna  E-field  is  not  aligned  with 
the  preferred  c-axis  alignment  is  not  as  apparent  (Fig. 
8b).  As  the  darkness  of  the  record  is  a  function  of 
the  voltage  amplitude  of  the  reflected  signal,  it  is  clear 
from  the  records  that  the  reflected  signal  strength  from 
the  ice  bottom  in  Figure  8b  was  weaker. 

Airborne  radar  sounding  measurements  were  also 
made  at  sites  Tig-79,  Boulder  Patch  (BP)  and  G  (Fig.  1). 
At  Tig-79  the  ice  was  covered  with  12  cm  of  snow. 

Ice  thickness  was  1 .59  m  and  the  local  c-axis  orienta¬ 
tion  was  150°  true.  The  direction  of  maximum  reflec¬ 
tion  amplitude  was  154°  true  (Fig.  9).  The  ice  exhib¬ 
ited  a  strong  coefficient  of  anisotropy  (4.7).  Not  far 
from  this  site,  at  the  1 978  field  season  location  Tigr3 
(Fig.  1 ),  the  under-ice  current  was  found  to  be  from 
130°  true  (Kovacs  and  Morey  1979a). 

The  ice  bottom  reflection  amplitude  data  from 
Site  BP  (shown  in  Fig.  1 0)  gave  a  coefficient  of  aniso¬ 
tropy  of  1.4.  The  direction  of  maximum  signal 
strength  was  1 35°  true,  in  agreement  with  the  preferred 
c-axis  orientation  of  the  bottom  ice,  which  was  found 
to  be  140°  true.  Snow  cover  at  the  site  was  less  than 
2  cm.  Ice  thickness  was  1.60  m.  Two  ice  cores  from 
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Figure  6.  Configuration  of  dual  antennas  mounted  to  the  helicopter  during  on-ice 
(a)  and  airborne  (b)  sounding  of  sea  Ice.  Vie w  a  also  shows  emergency  field  rotor 
deicing. 


Figure  7.  C-axis  direction  and  relative  voltage  amplitude 
of  the  wavelet  reflected  from  the  ice  bottom  vs  antenna 
E-field  azimuth  at  site  WD-79,  when  the  antenna  was  air¬ 
borne,  on  the  Ice  surface,  and  mounted  to  the  helicopter 
approximately  30  cm  above  the  ice  surface. 


b. 

Figure  8.  Graphic  records  of  airborne  impulse  radar  echo  sounding  of  sea  ice.  The  top  return  in  each 
record  Is  from  the  ke  surface  and  the  next  return  Is  from  the  ke  bottom.  The  long  period  variation 
In  the  record  is  due  to  changing  aircraft  altitude. 


Figure  9.  C-axis  direction  and  relative  volt-  Figure  10.  C-axis  direction  and  relative  volt¬ 
age  amplitude  of  the  wavelet  reflected  from  age  amplitude  of  the  wavelet  reflected  from 

the  ice  bottom  vs  antenna  E-field  azimuth  at  the  Ice  bottom  vs  antenna  E-field  azimuth  at 

site  Tig-79  (antehnd  airborne).  the  Boulder  Patch  site  (antenna  airborne). 


;  ?■ 


Figure  1 1.  Ice  cores  obtained  at  the  Boulder  Patch  site. 

this  site  are  shown  in  Figure  1 1 .  The  upper  three-  oriented  at  80°  true,  which  is  in  agreement  with  the 

quarters  of  the  ice  was  clear;  below  this  was  a  band  of  1 979  determination. 

dirty  ice  not  found  at  other  sites  studied  in  1 979,  and  At  site  G,  airborne  radar  soundings  were  made  from 

at  the  bottom  there  was  another  layer  of  clear  ice,  both  of  the  helicopter-mounted  antennas.  Each  anten- 

about  S  cm  thick.  na  was  operated  independently  in  the  transceive  mode. 

At  site  G,  the  snow  cover  was  less  than  1  cm  thick,  The  ice  bottom  reflection  amplitude  data  are  shown  in 

the  ice  was  1 .55  m  thick  and  the  preferred  c-axis  ori-  Figure  1 2.  The  difference  in  wavelet  amplitude  be- 

entation  was  80°  true.  In  1 977,  the  current  direction  tween  the  two  antennas  is  due,  in  part,  to  different 

at  this  site  was  measured  and  was  found  to  be  from  gain  settings.  Both  sets  of  antenna  data  reveal  that  the 

80°  true,  at  4  cm/s.  The  c-axis  was  also  found  to  be  ice  is  anisotropic,  with  the  maximum  reflection  occurring 


i  mf 
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in  the  same  direction,  i.e.  at  90°  true.  This  is  in  agree¬ 
ment  with  local  current  (Fig.  1 )  and  c-axis  direction 
observations.  The  latter  was  again  found  in  1979  to 
be  at  80°  true.  Also  shown  in  Figure  12  are  the  radar 
sounding  data  obtained  from  antenna  2  when  the  heli¬ 
copter  was  rotated  on  the  ice  surface.  These  data  are 
in  agreement  with  the  airborne  data  from  the  same  an¬ 
tenna.  The  fact  that  these  two  sets  of  data  have  simi¬ 
lar  amplitude  is  a  coincidence,  as  die  gain  setting  was 
lower  during  the  on-ice  measurements  than  during 
the  airborne  measurements. 

Our  measurements  at  the  freshwater  lake  site  were 
not  intended  to  detect  anisotropy  within  the  ice,  but 
were  made  for  control  purposes  and  as  a  check  to  vali¬ 
date  the  sea  ice  measurements.  However,  measure¬ 
ments  were  made  with  the  antenna  E-field  oriented 
at  0  and  90°  true.  The  resulting  radar  reflection  am¬ 
plitude  data  are  listed  in  Table  I. 

From  equations  previously  presented,  ice  surface 
and  bottom  coefficient  of  reflection,  apparent  surface 
and  bulk  dielectric  constants  and  bulk  velocity  calcula¬ 
tions  were  made  using  the  radar  reflection  amplitude 
and  signal  flight  times  vs  antenna  E-field  orientation. 
These  results  are  listed  in  Tables  I  and  II  for  the  1 978 
West  Dock  and  lake  ice  sites,  along  with  similar  data 
for  sea  ice  sites  1  and  3,  studied  in  1978.  The  latter 
site  information  is  listed  for  comparative  purposes 
and  because  the  bottom  coefficient  of  reflection  had 
not  previously  been  determined. 

From  the  elevated  antenna  data  the  average  surface 
reflection  coefficient  for  the  freshwater  lake  ice  (FW) 
was  determined  to  be  -0.277,  from  which  the  apparent 
dielectric  constant  is  calculated  to  be  3.12.  This  is  in 
agreement  with  accepted  published  values  (3.1  to  3.2) 
for  freshwater  ice.  The  effective  velocity  and  bulk  di¬ 
electric  constant  were  determined  to  be  0.170  m/ns 
and  3.1,  respectively.  Again,  these  values  are  in  agree¬ 
ment  with  accepted  published  values. 


Figure  12.  C-ax/s  direction  and  relative  voltage 
amplitude  of  the  wavelet  reflected  from  the  ice 
bottom  vs  antenna  E-fleld  azimuth  for  when  the 
antenna  or  antennas  were  "on"  the  Ice  and  air¬ 
borne  at  site  G.  From  the  polar  plots  It  appears 
that  antenna  2  has  more  directivity  then  anten¬ 
na  1. 


The  coefficient  of  reflection  from  the  lake  ice  bot¬ 
tom  was  determined  to  be  -0.206.  This  value  is  low  in 
comparison  with  the  theoretical  value  of  -0.684  (cal¬ 
culated  from  eq  6  using  3-1  for  e,  and  88  for  ew).  The 
lower  measured  value  for  the  coefficient  of  reflection 
is  believed  to  be  due  to  beum  spreading,  signal  absorp¬ 
tion  and  scattering  (A  and  5)  losses  within  the  lake  ice. 
The  A  and  5  losses  were  calculated  and  are  listed  in 
Table  III,  along  with  those  for  sea  ice  locations  WD-79 
and  sites  1  and  2.  The  loss  at  the  "sea  ice/water  inter¬ 
face"  was  calculated  using  12  for  e(  and  a  relative  com¬ 
plex  dielectric  constant  e*  of  310  for  ew  in  eq  6.  The 
relative  complex  dielectric  constant  of  the  seawater 
was  calculated  from: 

4- C,  VC  •*♦/,&-  « 

where  =  real  part  of  the  complex  dielectric  con¬ 
stant  of  seawater  ®  88  at  -1°C 
e"  =  imaginary  part  of  the  complex  dielectric 
constant 

o  -  conductivity  of  seawater  as  3  mhos/m  at 
-1°C 

f=  frequency  (175  MHz  from  Fig.  3a) 
e„  =  free  space  dielectric  constant  =  8.854x  10‘n 
farads/m. 

Therefore,  e*  =  88  +/308  and  |e*  |  *>  310,  which  is 
an  approximate  value  for  use  in  eq  6  to  take  into  ac¬ 
count  the  very  high  conductivity  of  seawater.  The  re¬ 
sulting  value  for  Rf/w  is  0.671 .  Because  the  electro¬ 
magnetic  boundary  at  the  sea  ice  bottom  has  a  higher 
brine  volume  (Fig.  1 3),  the  value  selected  for  e(  is  our 
best  estimate  at  this  time.  It  is  based  upon  the  extrapo¬ 
lation  of  the  curve  in  Figure  14  to  a  brine  volume  rep¬ 
resentative  of  that  found  at  the  bottom  of  sea  ice. 
Future  studies  are  planned  to  clarify  this. 
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(•)  Salinity  (%•) 

Figure  1.3.  Brine  volume  ami  salinity  distribution  vs  depth  at  site  WD-79. 


V 


Figure  14.  Average  effective  bulk  dielectric 
constant  vs  average  bulk  brine  volume  of  Ice. 


The  total  two-way  transmission  losses,  total  two- 
way  A  and  5  loss  and  the  one-way  A  and  5  loss  (from 
Table  III)  are  plotted  vs  the  average  bulk  brine  volume 
of  the  ice  for  each  site  (Table  IV)  in  Figure  1 S.  These 
losses  were  determined  from  the  elevated  antenna  data 
and  represent  the  minimum  loss  which  occurred  at  the 
maximum  recorded  voltage  amplitude  of  the  signal 
from  the  ice  bottom.  The  toss  in  the  sea  ice  is  shown 
to  increase  exponentially  with  increasing  average  bulk 
brine  volume.  It  is  apparent  that  when  the  average 
bulk  brine  volume  reaches  something  on  the  order  of 
UKf/m.A  and  S  losses  become  very  high  and  will 
probably  prevent  the  detection  of  the  sea  ice  bottom. 
This  would  apply  to  the  radar  system  used  when  it  is 
elevated  Urn  above  the  ice  surface  and  to  when  the 
antenna  E -field  is  aligned  with  the  preferred  c-axis  di¬ 
rection  of  the  ice.  Since  the  losses  become  much 
larger  when  the  antenna  E -field  is  oriented  perpend  icu- 


Table  IV.  Site  ice  thickness, 
average  bulk  brine  volume 
and  average  dielectric  con¬ 
stant. 


Site 

Ice 

thick. 

(mj 

Avg 

brine 

voi 

nui 

Avt 

diet. 

const. 

WD-79 

169 

47 

4.6 

1 

1-60 

76 

6.S 

3 

1  81 

66 

5.4 

FW 

201 

None 

31 

lar  to  the  c-axes  (because  the  voltage  amplitude  of  the 
reflected  wavelet  from  the  ice  bottom  decreases  signi¬ 
ficantly  in  this  direction),  it  is  even  more  unlikely  that 
the  sea  ice  bottom  could  be  detected  with  this  E-field 
alignment.  For  example,  at  this  E-field  direction  the 
total  two-way  losses  at  Sites  WD-79, 1  and  3  become 
-35,  -64  and  -60  db,  respectively.  These  values  are 
significantly  larger  than  the  losses  listed  in  Table  III. 
The  ice  bottom  might  still  be  detectable  if  the  antenna 
were  resting  on  the  surface,  as  the  spreading  losses 
would  be  less.  Conversely,  the  higher  the  antenna  is 
above  the  ice  surface,  the  lower  the  average  bulk  brine 
volume  of  the  ice  need  be  to  prevent  “seeing”  the  ice 
bottom.  Indeed,  in  the  spring  of  1978,  from  an  alti¬ 
tude  of  approximately  1 5  m,  the  ice  bottom  was  not 
detectable  at  site  1  (Fig.  1),  where  the  ice  had  an  aver¬ 
age  bulk  brine  volume  of  <s’76%»;  but  the  bottom  was 
intermittently  detectable  at  site  3  (Fig.  1),  where  the 
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Figure  15.  Transmission  tosses  in  ice  vs  its  average  bulk  brine  Figure  16.  Average  effective  coefficient  of  re¬ 
volume.  flection  from  the  ice  bottom  w  average  effec¬ 

tive  bulk  dielectric  constant  of  the  ice. 


average  brine  volume  of  the  ice  was  on  the  order  of 

66%o. 

The  average  effective  coefficient  of  reflection  at  the 
ice  bottom  for  when  the  antenna  is  elevated  vs  average 
bulk  dielectric  constant  (Table  I)  is  shown  in  Figure  16. 
The  average  effective  coefficient  of  reflection  is  diown 
to  become  smaller  as  the  average  effective  bulk  dielec¬ 
tric  constant  increases.  This  again  indicates  that  the 
ice  is  becoming  more  lossy  with  increasing  brine 
volume. 

Another  way  of  showing  this  is  by  the  plot  of  aver¬ 
age  effective  bulk  dielectric  constant  vs  average  bulk 
brine  volume,  as  shown  in  Figure  14.  In  this  graph,  e 
is  shown  to  increase  linearly  with  increasing  brine  vol¬ 
ume.  The  data  presented  are  from  the  1978  field  re¬ 


Figure  17.  A  verage  effective  velocity  of  the  electromag¬ 
netic  signal  in  ice  vs  Its  average  bulk  brine  volume.  The 
equation  for  the  line  passing  through  the  data  is  Ve  = 
0.17- 0.00068v  with  a  correlation  coefficient  of  0.90. 


suits  of  Kovacsand  Morey  (1979a)  and  from  the  1979 
field  season  for  ice  over  1 .6  m  thick.  Because  of  the 
limited  scatter  in  the  data,  the  second  parameter  should 
be  easily  obtainable  from  the  graph  after  one  parameter 
is  measured. 

The  data  in  Figure  14  can  be  replotted  as  the  average 
effective  velocity  vs  the  average  bulk  brine  volume,  using 
eq  2,  as  shown  in  Figure  17.  It  is  shown  that  the  aver¬ 
age  effective  velocity  decreases  linearly  as  the  average 
bulk  brine  volume  of  the  ice  increases. 

The  above  findings  are  significant,  for  they  provide  a 
way  to  not  only  measure  sea  ice  thickness,  but  also  to 
infer  its  strength  remotely  from  radio  echo  sounding  in¬ 
formation.  This  concept  has  also  been  suggested  by 
Rossiter  et  al.  (1977).  This  may  be  achieved  with  the 
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Figure  18.  Schematic  for  dual-antenna  electromag¬ 
netic  signal  flight  paths. 


use  of  a  dual-antenna  arrangement,  as  depicted  in  Fig¬ 
ure  18.  With  this  arrangement,  it  is  possible  to  deter¬ 
mine  the  effective  velocity  of  propagation  of  the  elec¬ 
tromagnetic  signal  in  die  sea  ice,  and  its  effective 
dielectric  constant,  irrespective  of  its  changing  proper¬ 
ties  due  to  temperature  or  thickness  variations.  One 
antenna  operates  in  a  transmit-receive  mode  and  the 
second  in  a  receive-only  mode.  The  antennas  are 
placed  a  fixed  distance  apart  on  the  ice  and  moved  as 
a  unit.  The  effective  propagation  velocity  ( Vt)  of  the 
signal  in  the  ice  is  then  determined  by: 


K 


(9) 


where  x  -  distance  between  centers  of  the  two  antennas 
rd  a  vertical  travel  time  from  transceiver  antenna 
to  and  from  subsurface  interface 
rx  =  travel  time  from  transceiver  antenna  to  sub¬ 
surface  interface  to  receive-only  antenna 
(rx)  plus  time  ra  =  air-wave  travel  time  be¬ 
tween  antennas  (Fig.  18). 


A  more  detailed  description  of  the  methodology  asso¬ 
ciated  with  use  of  a  dual-antenna  system,  and  some 
preliminary  results  which  show  that  ice  thicknesses 
measured  in  drill  holes  and  those  calculated  from  eq  9 
can  be  expected  to  be  within  10%,  can  be  found  in 
Kovacs  (1977)  and  Kovacs  and  Morey  (1979b). 

With  the  velocity  of  the  signal  thus  determined,  it 
is  posstole  to  determine  sea  ice  thickness  from  eq  2  and 
the  average  bulk  brine  volume  v  of  the  ice  by: 


0-17-  Ve 
0.0007 


(10) 


where  eq  10  fits  die  line  passing  through  the  data  in 
Figure  17. 

Since  the  mechanical  properties  of  sea  ice  are  a  func¬ 
tion  of  its  brine  volume,  it  follows  that  a  strength  prop¬ 
erty  (o)  of  the  ice  may  now  be  estimated  remotely  by 
radio  echo  sounding.  Most  strength  properties  of  sea 


ice  can  be  represented  by  an  equation  of  the  form: 

o=a-b^v  (11) 

where  a  and  b  are  constants.  For  example,  Vaudrey 
(1977)  found  the  following  relationship  from  beam 
test  results: 

of  =  9.8- 0.62  v/P  (12) 

where  of  is  the  flexural  strength  of  sea  ice  in  kg/cm2 . 
Thus,  with  the  apparent  thickness  and  strength  of  sea 
ice  known,  it  follows  that  the  bearing  capacity  of  the 
ice  may  also  be  approximated  by  analytical  methods 
which  are  beyond  the  scope  of  this  report,  e.g.  see 
Vaudrey  and  Katona  (1974),  Frederking  and  Gold 
(1976),  Nevel  (1978, 1979)  and  Johnson  (1980). 


CONCLUSIONS 

This  study  supports  previous  work  (Kovacs  and 
Morey  1978, 1979a)  which  revealed  that  for  sea  ice 
with  a  bottom  structure  in  which  the  horizontal  c-axes 
of  the  ice  crystal  platelets  are  aligned,  this  ordered  bot¬ 
tom  structure  is  an  effective  polarizer  of  transverse 
electromagnetic  waves.  In  short,  when  a  linearly  polar¬ 
ized  antenna  E-field  is  aligned  parallel  with  the  pre¬ 
ferred  c-axis  azimuth,  a  maximum  signal  return  is  re¬ 
ceived  from  the  ice  bottom,  but  when  the  antenna  is 
aligned  perpendicular  to  the  preferred  c-axis  orienta¬ 
tion,  the  signal  is  significantly  reduced  or  eliminated. 
This  effect  has  been  attributed  to  the  ordered  arrange¬ 
ment  of  the  brine  inclusions,  which  are  believed  to 
create  a  unique  array  of  lossy,  parallel  plate  waveguides 
at  the  ice  bottom  (Kovacs  and  Morey  1978).  Recent 
laboratory  studies  of  artificial  dielectrics  with  polariz¬ 
ing  effects  similar  to  sea  ice  support  the  concept  of  a 
sea  ice  model  in  which  the  ice  bottom  is  composed  of 
such  an  array  of  waveguides  (Morey  and  Kovacs,  in 
prep.). 

In  this  study  additional  data  are  presented  which 
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also  verify  that  the  electromagnetic -dependent  proper¬ 
ties  of  sea  ice  vary  in  the  horizontal  plane,  as  does  the 
anisotropy.  It  is  now  known  that: 

1 .  The  average  effective  bulk  dielectric  constant, 
and  therefore  the  average  effective  velocity  of  the 
electromagnetic  pulse  used  in  this  study,  is  dependent 
upon  the  average  bulk  brine  volume  of  the  sea  ice. 

2.  Sea  ice  anisotropy  may  be  detected  by  radio  echo 
sounding  measurements  made  not  only  on  the  ice  sur¬ 
face  but  also  from  an  airborne  platform. 

3.  The  effective  coefficient  of  reflection  from  the 
sea  ice  bottom: 

a.  decreases  exponentially  with  increasing  aver¬ 
age  effective  bulk  dielectric  constant  of  the  ice, 

b.  decreases  with  increasing  bulk  brine  volume, 
and 

c.  is  typically  one  to  two  orders  of  magnitude 
lower  than  the  coefficient  of  reflection  from  the 
ice  surface. 

4.  The  losses  in  sea  ice  increase  “exponentially” 
with  increasing  average  bulk  brine  volume. 

5.  Snow  on  the  ice  reduces  the  voltage  amplitude 
of  the  electromagnetic  wavelet  reflected  from  the  ice 
bottom. 

The  fundamental  relationship  between  the  average 
bulk  brine  volume  of  sea  ice  and  its  electrical  and 
strength  properties  is  also  discussed.  It  is  shown  that, 
in  principle,  it  should  be  possible  to  not  only  determine 
sea  ice  thickness,  but  also  to  estimate  its  strength  re¬ 
motely  with  the  use  of  a  properly  designed  electromag¬ 
netic  sounding  system. 

Implementation  of  the  dual-antenna  mode  for  de¬ 
termining  velocity  of  electromagnetic  wave  propaga¬ 
tion  in  sea  ice  may  be  limited,  especially  in  an  air¬ 
borne  mode,  to  flight  altitudes  of  less  than  twice  the 
antenna  separation,  because  of  inherent  limitations  in 
accurately  measuring  signal  flight  times  of  less  than 
about  0.5  ns.  Also,  this  method  assumes  that  most  of 
the  electromagnetic  anisotropic  effects  are  occurring 
in  the  lower  10%  of  the  ice  sheet,  as  postulated  in 
Kovacs  and  Morey  (1979a)  and  Golden  and  Ackley 
(in  press).  For  example,  eq  9  is  based  upon  straight 
ray  path  propagation  in  a  homogeneous  dielectric, 
which  sea  ice  is  not.  The  propagation  distance  through 
sea  ice  will  not  be  a  straight  line  but  will  be  bent, 
based  upon  Snell's  law  for  a  material  with  increasing 
refractive  index  (which  appears  significant  for  growing 
winter  sea  ice  only  near  the  ice  bottom).  Preliminary 
results  with  the  use  of  the  dual  antenna  system  indi¬ 
cate  that  ray  path  bending  is  not  a  serious  problem 
(Kovacs  1977).  Further  evaluation  of  this  is  planned. 
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APPENDIX:  DATA  ANALYSIS  PROCEDURES 


The  radar  data  were  recorded  on  an  analog  magnetic 
tape  recorder  after  down-conversion  in  the  radar  re¬ 
ceiver.  The  down-conversion  consists  of  a  sampling 
process,  wherein  the  nanosecond  (ns)  time  frame  is 
converted  to  a  millisecond  (ms)  time  frame  for  record¬ 
ing  and  display.  A  typical  scan  iishown  in  Figure  A! . 
The  voltage-amplitude  vs  time  plot  consists  of  a  series 
of  wavelets  representing  reflections  from  the  indicated 
interfaces,  e.g.  the  air/ice  interface.  The  radar  receiver 
contains  a  time-variable-gain  (TVG)  circuit  which  am¬ 
plifies  the  signal  as  time  increases.  The  slope  of  the 
TVG  function  was  also  recorded  on  tape.  During  data 


analysis,  as  when  the  relative  peak -to -peak  voltage 
levels  of  various  reflected  wavelets  were  being  com¬ 
pared,  the  effect  of  the  TVG  function  was  removed. 
Also,  part  of  the  field  calibration  procedure  requires 
calibrating  the  time  base  of  the  radar,  using,  for  ex¬ 
ample,  a  length  of  coaxial  cable  exactly  10  ns  long,  or 
by  suspending  the  antenna  a  measured  distance  above 
a  large  metal  reflector. 

A  digital  signal  analyzer  was  used  for  laboratory 
evaluation  of  the  radar  data.  The  same  tape  recorder 
used  in  the  field  to  record  the  data  was  used  to  play 
the  data  back  through  the  analyzer.  An  average  of  ten 
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Figure  Al.  Representative  scons  obtained  from  sea  fee  sounding  using  Impulse 
radar  vdten  antenna  was  on  the  Ice  surface  (top)  and  elevated  above  the  ke  sur¬ 
face  (bottom). 
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Figure  A2.  Representative  reflected  wavelet  frequency  spectrum  as  generated 
by  digital  signal  analyzer. 


scans  was  created  to  reduce  the  effect  of  random 
noise.  Movable  cursors  in  the  analyzer  were  positioned 
at  the  same  point  on  each  wavelet,  e.g.  the  first  zero 
crossing  of  the  air/ice  and  ice/water  interface  reflec¬ 
tion,  and  the  time  difference  in  milliseconds  recorded. 

A  conversion  factor  was  calculated  using  the  10-ns 
period  between  the  calibration  pulses  as  displayed  on 
the  digital  signal  analyzer.  For  this  example  the  pulse 
period  on  the  analyzer  is  33  milliseconds;  therefore 
the  conversion  factor  is  3.3x10'* 4-1  O'*,  or  3.3x10*. 

The  conversion  factor  was  used  to  calculate  the  two-way 
travel  time  in  nanoseconds  between  two  wavelets  as 
displayed  on  the  analyzer. 


The  digital  signal  analyzer  can  also  calculate  t>*  fre¬ 
quency  spectrum  of  a  wavelet.  This  was  done  by  digi¬ 
tizing  only  the  reflected  wayelet  from  a  single  interface 
and  then  calculating  the  frequency  spectrum  of  that 
wavelet.  The  spectrum  is  displayed  as  the  log  amplitude 
(in  db)  vs  frequency  (in  Hz),  as  shown  in  Figure  A2. 

The  -3db  points  of  the  spectrum  were  determined  and 
the  center  frequency  calculated  using 

'c’VhOI*- 

For  example,  fc  in  Figure  A2  is  67  Hz.  This  value  is 
then  multiplied  by  the  conversion  factor  to  determine 
the  actual  center  frequency,  i.e.  3.3x1 0‘x67  Hz=220 
MHz. 
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